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^ . Abstract 

TT*- The contact process is a stochastic process which exhibits a continuous, absorbing-state 

ry^ I phase transition in the Directed Percolation (DP) universality class. In this work, we con- 

sider a contact process with a bias in conjunction with an active wall. This model exhibits 
waves of activity emanating from the active wall and, when the system is supercritical, 
propagating indefinitely as travelling (Fisher) waves. In the subcritical phase the activity 
1-^ I is localised near the wall. We study the phase transition numerically and show that certain 

properties of the system, notably the wave velocity, are discontinuous across the transition. 
Using a modified Fisher equation to model the system we elucidate the mechanism by which 
the the discontinuity arises. Furthermore we establish relations between properties of the 
travelling wave and DP critical exponents. 
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^ ! 1 Introduction 

The contact process [1] exhibits a nonequilbrium phase transition from an absorbing state, 

'^ I where the system ends up in an inactive configuration, to an active, fluctuating state. Although 

Q ■ originally introduced as a microscopic model for epidemic spreading, this lattice model and 

O . its relatives have been used to describe a variety of systems including percolation, wetting, 

reaction-diffusion processes, branching and annihilating random walks [2] and phase transitions 

CN ■ in more exotic settings, such as between two turbulent states in nematic liquid crystals [3], 

and in proliferating microbial populations under gravity [1]. The continuous phase transition 

[~^ I from the inactive to active phases falls into the Directed Percolation (DP) universality class 

r^ ' which is thought to pertain for any such microscopic model exhibiting an absorbing state phase 

transition, in the absence of any additional symmetries or conservation laws [5l|6]. 

The nature of the transition may be modified by the introduction of an active boundary 

t^ I into the system that ensures that the activity in the system never dies out [TlIlOj. These 

studies have mainly focused on surface critical behaviour and the emergence of new surface 

critical exponents. In the counterpart to the inactive phase the activity is confined to the 

K^ \ boundary region whereas in the counterpart to the active phase the activity spreads from the 

^ • boundary through the system. The spreading of activity through the system from the boundary 
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is conveniently illustrated in the mean-field description of the DP universality class [2] which has 
the same form as the Fisher-KPFU equation [IT]. This nonlinear partial differential equation 
exhibits travelling wave solutions in which the active phase invades inactive regions with a 
well-defined velocity for the domain wall that separates the two regions. 

In this work we consider a variation of the contact process with the two features of an active 
boundary and advection away from the boundary which we refer to as the Driven Asymmetric 



^Here Fisher-KPP stands for Fisher-Kolmogorov-Petrovskii-Piskunov. It is alternatively called the Fisher- 
Kolmogorov or Fisher equation. 



Contact Process (DACP). The Asymmetric Contact Process has been previously studied in the 
absence of boundaries and yields a phase transition in the DP universality class, as expected 
|12H14j . The addition of an active boundary, which serves to drive the system, then modifies 
the nature of the transition as described above [9]- However, what we discover is that the 
DP transition now has a discontinuous aspect in the sense that the velocity of the wavefront 
which carries the activity from the boundary jumps discontinuously at the critical point. Thus, 
intriguingly, the continuous DP transition is accompanied by a velocity discontinuity in the 
presence of an active boundary and advection. 

There has been recent interest in the addition of an advection term to the Fisher-KPP 
equation in the presence of a boundary since a Galilean transformation no longer serves to 
remove the advection term (as would be the case in the absence of boundaries). When the 
advection is directed towards the boundary there is a competition between the advective velocity 
and the Fisher wave velocity, leading to a phase transition [^ 1151 116] from a phase with activity 
localised near the boundary to one where the Fisher wave invades the whole system. The 
addition of noise into this system leads to a more complicated scenario where, in the case of 
a reflecting boundary studied in [TTj, the low activity phase could either be localised near the 
boundary or be the absorbing state. However, the discontinuous transition observed in that 
work where the bulk density jumps at the transition, is distinct from that studied in the present 
work. 

To understand the transition we observe, we consider a phenomenological mean-field de- 
scription of the system in the form of a Fisher-KPP equation which includes the effects of 
asymmetry in the form of an advective term. We show that the solution can be thought of as 
a Fisher wave moving within an envelope given by the stationary density profile. Thus, in the 
subcritical regime the system exhibits attenuated waves with non-zero velocity but whose am- 
plitude decays to zero away from the boundary. In the supercritical regime, on the other hand, 
the stationary profile has a non-zero limit far from the boundary, thus the Fisher wave propa- 
gates into the bulk with a constant amplitude and non-zero velocity. This is the mechanism for 
the observed discontinuity of the velocity. 

The paper is organised as follows: In section 2 we define the microscopic model and in 
section 3 we present numerical evidence for the discontinuous phase transition. In section 4 we 
discuss the observed scaling behaviour and identify the Directed Percolation critical exponents. 
In section 5 we discuss a mean field theory in the form of a modified Fisher-KPP equation. 
Using the mean field picture we then re-examine in section 6 the simulation results to show that 
a similar picture holds. We conclude in section 7. 

2 Microscopic model 

The driven asymmetric contact process (hereafter, DACP) is a stochastic model defined on 
a one-dimensional lattice. A microscopic configuration is specified by the set of occupation 
numbers {T(i)} where Tt{i) = (1) indicates that site i is inactive (active) at time t. The 
leftmost site of the system is kept permanently active, Tt(0) = 1 Vt, making the process driven. 
Active sites can activate inactive sites directly to the right of them at rate r or become inactive 
with rate 1. See Fig. [T] for an illustration. 

The driving from the left boundary eliminates the absorbing state associated with the basic 
contact process. However as the rate r goes through a critical value re there is still a continuous 
phase transition from a phase where the active sites can only spread a finite distance from the 
left boundary to a phase where they can spread to infinity. See Fig. [2] for an illustration of the 
density profile in the two states. 





Figure 1: The Driven Asymmetric Contact Process. In this model, the first site is kept active. 
Active sites can activate inactive sites to the right of them at rate r or become inactive 
with rate 1. 




Figure 2: A sketch of the density profile of the two states. In the subcritical state (red) the 
density decays exponentially whilst in the supercritical it propagates at a constant 
bulk density. 



An interesting and attractive feature of the model is that the total asymmetry of the acti- 
vation dynamics means that the state of site A'^i is independent from that of N2 for N2 > Ni. 
This implies that any calculation or simulation performed on a finite system of A^ sites gives 
the exact behaviour for the first N sites of an infinite system. In other words, the introduction 
of a right boundary does not introduce any finite-size effects. 

We now summarise previous work on asymmetric contact processes. Three mathematical 
papers |12fll4] have looked at how introducing asymmetry into the contact process (without 
a boundary drive) leads to the emergence of a second order parameter, the probability that 
the origin is active as t —t- 00, in addition to the probability that the process remains active 
indefinitely. At total symmetry the two order parameters coincide and at total asymmetry the 
second disappears. An approximate analytical and numerical study of the totally asymmetric 
contact process (without a boundary drive) [18] found a continuous transition at the critical 
rate Tc = 3.306(4) with the critical exponent f3 = 0.2760(1) which is in agreement with the value 
for the DP order parameter exponent in one spatial dimension. In addition a two-dimensional 
generalisation of the DACP was constructed [9j and used to study wetting and interface phe- 
nomena 1191. 



3 Numerical evidence for a discontinuous velocity transition 

We begin our study of the DACP by presenting numerical data that suggests the velocity of 
the wave emanating from the left boundary is discontinuous at the phase transition point. 



3.1 Simulation details 

We performed direct Monte Carlo simulations of the microscopic DACP dynamics specified in 
section 2 above. Each run was initialised with the lattice unoccupied, apart from the site 
which is kept permanently active. During the simulation a list of active sites is maintained. In 




Activation Rate (r) 

Figure 3: A plot of the velocity vs. activation rate for different run times. The two shorter 
runs spanned activation rates in the interval re [1,8] and the longer two re [1, 4], in 
both cases the intervals were subdivided into 500 equally spaced values of r for which 
the simulation was run. Note that as the simulation time is increased a discontinuity 
appears to develop. 



each elementary update, one of the n active sites is chosen from this list at random. With prob- 
ability ^7^, its right neighbour is activated; otherwise the chosen active site becomes inactive. 
Since our aim is to measure a velocity, it is important to keep track of the length of time 
associated with each update. Since we attempt an update on any of the n active sites the total 
attempt rate is A = n(l + r). In principle, the size of the time increment. At, for each step 
should be sampled from an exponential distribution with mean 1/A. However since the time 
increments of any realisation are small the distribution becomes sharply peaked and we can 
instead use the mean of the distribution as the size of our timestep. This approximation is 
computationally favourable and we have verified that it made no difference to the results. Thus 

, where n is the total number of active 
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n(l+r) 



we take the time taken for each update to be At 
sites before the update takes place. 

The front velocity of a single run was defined as the position of the rightmost active site at 
the end of the run, divided by the time taken for the run to end. This definition was chosen 
since it is unambiguous in the microscopic model. The simulations were run until a prespecified 
amount of time had elapsed: the lattice expanded as needed to accommodate all active sites. 
To reduce the noise in the data the simulation was repeated for 500 runs for each value of r and 
the quantities of interested averaged over this ensemble of runs. 



3.2 Results 

The first interesting result pertains to the asymptotic front velocity. The simulation was run 
for several different simulation times as shown in Fig. [3l Our definition of the front velocity 
implicitly assumes that an asymptotically constant velocity is reached on a timescale short 
compared to that at which the simulation ends. We anticipate that this may not be the case for 
values of r near the critical value re, since relaxation times diverge here. Therefore measurements 
of the velocity near Tc are subject to finite-time corrections. As can be seen from Fig. [3l 
increased run time shows a sharpening of the velocity as a function of r near re, suggestive of 
a discontinuous transition in the infinite-time limit. 

We also studied the density profile in the two different phases. The position and shape of the 
front can be investigated by looking at the density profile for a single activation rate at different 
times. Choosing two values of r close to, and on either side of, the critical rate we observe two 
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Figure 4: (a) Density vs. position in the supercritical regime at different times. A wave of 
constant bulk density invades the empty lattice as time progresses, (b) Correspond- 
ing data for the subcritical regime. The wave is attenuated by an envelope fitted 
empirically by the solid line. 

very different scenarios. For the supercritical case (Fig. H^) there is a wave of constant bulk 
density invading the empty lattice. The front simply propagates away from the boundary at a 
constant velocity. 

For the subcritical case (Fig. (Da) we observe an attenuated wave: the front propagates away 
from the boundary whilst at the same time decaying away. However the decaying density seems 
to follow an envelope of sorts. As a guide to the eye this envelope has been added to the graph. 
We will return to the idea of this envelope again in section [SJ 

It is perhaps surprising at a first glance that the continuous DP transition should be ac- 
companied by a discontinuity in the front velocity. In the remainder of this work, we elucidate 
the mechanism behind this discontinuity. 

4 Directed Percolation scaling picture 

Despite the absence of an absorbing state in the DACP, we show in this section that the 
behaviour of the DACP described in the previous section can in fact be interpreted within the 
universal scaling picture associated with the Directed Percolation phase transition. 



4.1 Steady-state density profile 

We shall consider the steady state density profile of the system pj = lirat^ooiTj) ■ A suitable 
order parameter for the system is the steady-state density as j — )• c«, which we denote by 'p. In 
the supercritical phase, close to criticality, we expect the order parameter to scale as p ~ A'^, 
where A = |r — rd is the distance from criticality and f3 is the order parameter exponent which 
we expect to be equal to the value for DP in one spatial dimension. We can also define a 
characteristic length scale ^ as the steady-state density decay length in the subcritical phase 
where pj ~ e~^'^. Then, since the length scale diverges at criticality we expect it to scale as 
S, ~ A~'^ near criticality where i/ is a correlation length exponent. At criticality we expect the 
profile to decay as a power law with exponent 6. 

Near criticality we further expect the scaling form pj ~ j~ g{j H) with g{u) a scaling 
function obeying \mi^^^ g{j / £^) = g{0) = constant such that at criticality pj ~ j~ . However 
since the steady-state density approaches a non-zero constant as x — )• oo in the supercritical 
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Figure 5: Collapse of the steady-state density profile pj for a range of activation rates r onto two 
master curves, one for the supercritical regime (upper curve) and one for the subcrit- 
ical regime (lower curve). The collapse was obtained by varying the free parameters 
Tc, i5 and v. Figure reproduced from 0. 

In [9], the critical point r^ and exponents 5 and v were obtained by plotting pjj against u = 
j/^ = j\r — r^ for different values of r, and varying r^^ b and v until the best data collapse (as 
judged by eye) was obtained. The resulting collapse, obtained for r^ = 3.3055(5), 5 = 0.1640(5) 
and V = 1.7(2), is reproduced from [9] in Fig. [5l 

We compare these measurements with the established values of the DP exponents in table [TJ 
We observe good agreement with the DP values that apply in one spatial dimension (ID) as 
long as we identify the exponent v with the DP temporal correlation exponent vh, as opposed 
to the distinct, and independent, spatial exponent iy±: we discuss this point in more detail 
shortly. For future reference, we have also included the exponents obtained within a mean-field 
approximation in table [TJ For an in-depth study of DP exponents and scaling we refer the 
reader to [2]. Here, we focus more on the dynamic behaviour of the DACP and the behaviour 
of the active front as it moves out from the left boundary. 

DP (MF) [2] DP (ID) [20J DACP [9] 



5 



1 
1 

1/2 



0.159464(6) 
1.733847(6) 
1.096854(4) 



0.1640(5) 
1.7(2) 



Table 1: Critical exponents for Directed Percolation (DP) in a mean field approximation (MF) 
and one spatial dimension (ID), and for the Driven Asymmetric Contact Process 
(DACP) of the present work. 



4.2 Advection dynamics: the shearing of DACP into DP 

The easiest way to understand the effect of asymmetry (or advection) in the contact process is to 
directly compare a space-time plot of the DACP dynamics with its contact process counterpart 
and its description in terms of DP scaling exponents — see Fig. [6j From these plots we identify 
two characteristic angles. In the DACP, Fig. [6^, the activity emanates from the wall with an 
axis that is an angle to the vertical. This axis corresponds to the time direction in the basic 
contact process. Fig. [6)3. The leading edge of the activity emerges at an angle (p to the DP time 
axis. We thus picture the DACP as a spatial shearing of DP by the shear element tan^. One 
consequence of this shearing is that the spatial correlation length measured in the steady-state 
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Figure 6: Space time plots of (a) the DACP and (b) DP. In DP, a cone of activity spreads at 
an angle to the time (vertical) axis; to obtain the DACP, this cone is sheared by a 
further angle 6 due to the advection. Note that the activation rates have been chosen 
to allow easy identification of the two angles in these figures: the distance from the 
critical point is different in the two cases, and hence there is no significance of the 
density of the DP clusters being greater than that for the DACP. 



density profile of the DACP is a linear combination of the temporal and spatial correlations 
lengths in DP. Such that Cj_ = C± + ClI tan 9 and ^f, = ^ii , with primed variables pertaining to 
the DACP. Hence, we expect the DACP correlation length to diverge with the faster-growing 



DP correlation length as r 



and therefore that u should be identified with the larger of 



the DP exponents (i.e., i^y, as above). 

We now turn to the dynamics. Above the critical point in DP, an activity wave travels at a 
velocity udp = tsii(j). If, near criticality, there is a single characteristic length and time scale, 
this velocity must be given by their ratio. That is 



vjiv = tan ' 






A-^ with X 



11 - ^^ 



(1) 



Since v\\ > u±, this velocity vanishes as r —t- re from above. Below criticality, activity does not 
spread out indefinitely, and so the wave velocity is zero when r < Vc- Thus, the wave velocity 
is continuous across the transition in the absence of advection. 



We obtain the supercritical wave velocity for the DACP, 
That is, 

^^ sin (j) + tan cos 



"sup 



, by applying the shear tan^. 



i^sup 



COS( 



tan (j) + tan 9 = udp + s = s + j4A-^ 



(2) 



where we have introduced the advection s = tan 6. (There may be subleading corrections coming 
from corrections to scaling in the relation ([1]).) The actual value of s is a nontrivial emergent 
property of the stochastic DACP dynamics. However, a key point is that it can take a nonzero 
value in the subcritical regime which yields the intrinsic velocity VsuW- activity still propagates 
at rate Usub = s from the active boundary, despite the fact that it dies out after some finite 
time. Thus, the apparent (i.e. the observed) wave velocity, measured far from the origin, may 
jump discontinuously from zero to some nonzero value Vc = s{rc) at the critical point. This 
is in contrast to the DP wave velocity which goes to zero at criticality, as shown above. The 
leading-order behaviour of the supercritical wave velocity in the DACP is then v. 



sup 



Vc + ^A^ 

where X = ^11 ~ ^±- We remark that this expression provides a means to measure the smaller 
correlation length exponent, z/_l, which cannot be accessed from the steady-state density profile 
alone. 



5 Mean-field theory 

Having proposed a mechanism for the discontinuous velocity transition observed in section 3, 
we now demonstrate exphcitly that it is at work within the mean-field formulation of the DACP 
dynamics. 

We know that for a DP model, such as the contact process, the continuous-space mean-field 
description is given by a Fisher-KPP equation. The Fisher-KPP equation was first introduced 
in [2H - I23J . see |llj for a review. We begin with the general form of the Fisher-KPP equation 
supplemented with an advection term to take into account the boundary drive present in the 
DACP: 

dp o 2 ,dp , r,^'^P /Q^ 

At this stage a, [3, D and r' are phenomenological parameters. We will match these up with the 
parameters of the microscopic model in section 16.11 below. Meanwhile we note that we expect 
r' and D to be smooth functions of the microscopic advection coefficient r and that a (x r — Vc- 
We first establish the two distinct stationary regimes which are determined by the sign of 
a. If we set -g^ = in ([SD we find the stationary density p*{x) is given by 

n * o *2 /dp* „9 p* 

= ap -Pp -r^ + D^. (4) 

The solution takes the following forms for large x. 

For a < {subcritical case) : a density profile which decays exponentially to zero and has 
the form for large x, p^^^^ix) — Ae~ 



For a > {supercritical case) : a density profile which decays exponentially to a non-zero 



value and has the form for large x, Psup{x) — p + Be "^^ with p 



In both cases the decay constant may be written as 

-r' + (r'^ + AlalDV^^ 
^ = 2D ■ (^) 

For small |a|, A = 0(|a|). Thus the characteristic lengthscale ^ diverges as |q;|~^ at criticality. 
Following the arguments of section 4, ^ is expected to diverge with the DP exponent z^ii. From 
table [H we see indeed that uu = 1 in the mean field. 

Based on the envelope observed in the simulation of the stochastic system (see section 3) we 
assume that the full time-dependent density can be described by p{x,t) = p*{x)f{x,t) where 
p*{x) is the stationary solution to (jl]). Equation 1^ then becomes 

p*f = ap*f - Pp*^f - r'{p*'f + p*f') + D{p*"f + 2/7' + P*!") ■ (6) 

Dividing through by p* and using the definition of /?*, dU, to eliminate a we obtain a modified 
Fisher-KPP equation for the wave /(x, t) that sits inside the envelope: 

df df d'^f 

In this equation, the advective velocity of the modified wave, s, is given by 

s = r —- . (8) 



Note that this advective velocity, s, is generahy increased over the bare quantity r' as a con- 
sequence of the envelope. The coefficient, f3p*{x) of the non linear growth term /(I — /) is 
generally x dependent, but in the supercritical phase it decreases to a constant value for large 
X thus recovering the Fisher-KPP equation far from the boundary. On the other hand, in the 
subcritical regime the non-linear term decreases to zero far from the boundary. 

First we consider the supercritical phase. Far from the boundary ([7]) becomes, using p* — t- 
a/f3 and s ^ r' 

which is the usual Fisher wave equation with advective coefficient r'. Following the usual 
approach (see e.g. [TT]), one assumes a travelling wave form f{x,t) = f{x — vt) = f{z) and 
linearises for large z where / is small. The solution is of the form / = e~^^ with 

iv - r') ± J(v - r'Y - ADa , , 



Since /i € M_|_ we require {v — r')^ > ADa ^ v > r' + \^ADa = r^min- We now assume (as for 
the normal Fisher wave starting from an initial sharp front) that v = fmin- Thus 



i;,up = / + VId^ . (11) 

Of course ([TT]) could simply be obtained by adding the advective velocity r' to the usual Fisher 
wave velocity in the absence of advection, (AaD) ' . 

In the subcritical phase, the modified wave described by ([7]) has a velocity given by a distinct 
expression. Here, for large x, ([7]) becomes 

/ = A/3e-^V(l - /) - sf + Df" (12) 

where 

s = r' + 2DX = {r'^ + A\a\Dy/'^ . (13) 

Note that the advective velocity of the attenuated wave is increased over the value r'. 

The spatial dependence of the coefficient of the non-linear term makes the analysis of equa- 
tion p2|) non-trivial. Here we content ourselves with a heuristic picture. Initially the presence 
of the nonlinear term will mean that a nonlinear travelling wave emanates from the boundary. 
However as the front moves away from the boundary the nonlinear term becomes less important 
and we expect the front to broaden and the velocity to decrease. Finally as the front of the 
wave moves further away and x ^ 1/A the equation for / reduces to a diffusion equation with 
advection 

/ = -sf + Df" . (14) 

The wavefront thus broadens diffusively over time. Therefore, at late times, the modified density 
profile / will be approximately that given by a diffusion equation. This diffusive front moves 
with velocity 

Vsnh = s = {r'^ + 4\a\D)'/\ (15) 



and the profile itself takes the form 



f{x,t) ~ -erfc 



X — st 



2VDt 



(16) 



Thus the width of the front is ultimately \'Dt. 



To summarise, the analysis of the phenomenological Fisher-KPP equation ([3]) shows that 
there are two possible regimes, according to the sign of a. When a is negative, the density 
decays to zero as it moves away from the origin; when positive, it propagates away from the 
origin with a constant bulk density at a constant velocity. The dynamics can be couched in 
terms of a modified wave, fix — vt), travelling within the envelope of the steady-state density 
profile, p*{x). This wave is governed by the modified Fisher-KPP equation, ([7|), and has a 
velocity in these two states of t'sub = Vr''^ — 4:Da and iigup = r' -\- \J \.Da. Note that v is 
continuous at the transition a = although its derivative is not. 

In the subcritical regime, this velocity coincides with the advective velocity s introduced in 
section 4 in terms of the angle of shear from DP to the DACP. As noted previously, one effect 
of the envelope is to force the apparent velocity (i.e., that observed in simulations) to decay to 
zero: activity can probe only a finite distance from the origin as time t — t- oo. Above the critical 
point, the observed velocity is Wsup and so, across the transition, this observed velocity exhibits 
a discontinuity. We note from the form of t^sup; the velocity grows from its critical value as 
■u — f c ~ a-^ where x = 1/2. This value of x agrees with the prediction x = v\\— i^± from the 
mean-field theory for DP (see table [TJ . 

6 Modified travelling wave in the stochastic DACP 

We now revisit the simulations of the DACP in the light of what we have learnt from the mean- 
field theory. We first express the phenomenological parameters appearing in the Fisher-KPP 
equation, Q, in terms of the stochastic activation rate r. We then examine more closely the 
numerical evidence for the picture of an modified wave travelling within the stationary density 
profile, and in particular, the prediction for the growth of the wave velocity just above criticality, 
i.e., V — Vc ^ a^ where x = z^y — '^± given by the appropriate DP exponents. 

6.1 Identification of the mean-field phenomenological parameters 

We first revisit the mean-field theory of Section [5] and compare to the microscopic DACP dy- 
namics. One may make and make a heuristic identification of the phenomenological parameters 
in ([3]) by considering first of all an exact equation for the rate of change of density at site i in 
the DACP. This reads 

f^init)) = r{r,^,{t) [1-rM) " {n{t)) , (17) 

where the first term comes from site i — 1 activating site i at rate r if the former is active and 
the latter inactive, and the second term from the decay of site i at unit rate when it is active. 
The mean-field approximation is to write {Ti{t)Tj{t)) = pi{t)pj{t), where pi{t) = (tj). 

Pi = rpi^i{l- Pi)- Pi Vi > 1 
Po = l Vt 

The steady state solution of this spatially discrete equation is provided in the appendix. 

Here we move over to continuous space and expanding to second order spatial derivatives, 
we find, a = r — Tc, /3 = r, r' = r(l — p) and D = ^, where Tc = 1. This suggests a density 
dependence in r' and D: however, as we now explain, this is not expected to affect the wavefront 
behaviour. 

In the subcritical regime the density profile tends to zero far from the boundary therefore 
r' — 7- r and D — t- r/2. Thus 

^sub = yf'^ — 2ra = y/r'^ — o? ~ rg — V ■ ■ ■ (19) 
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Figure 7: (a) The density vs. position with the envelope divided out (r = 3.3). (b) The density 
vs. position/time. The collapse shows that the velocity remains constant. 

In the supercritical regime, on the other hand, p will tend to zero at the tip of the wavefront 
but be nonzero behind the Fisher wavefront. Thus at the tip r' = r and the supercritical Fisher 
wave velocity becomes 

+ \/2ra — Tc + \j2rca . 



•^sup 



(20) 



Therefore the derivative of the modified wave velocity is discontinuous at the phase transition. 
It is also interesting to note that with our identification of r' and /3, behind the wavefront in 
the supercritical regime we have r' = r(l — a/r) = r^ and the effective advection is fixed at the 
critical value. 



6.2 Stochastic modified wave dynamics 

Returning now to the stochastic simulations of the DACP, we investigate first the picture 
of a modified wave travelling within a density envelope. In Fig. IDa, the density profile in 
the subcritical phase is shown at different time points along with an envelope of the form 
p*{x) = exp(— ^02; — Aix '^), with the fitting parameters Ai all positive. As we get closer to 
criticality A2 decreases, indicating an approach to the scaling form e~^'^x~ . 

To obtain the modified wave f{x — vt), we divide the numerical density profiles by this 
envelope equation. The result of this procedure is shown in Fig. [7^, which clearly shows a wave 
with constant bulk density invading an empty lattice. The constancy of the wave velocity can 
be checked by dividing the x coordinate by time (Fig. U}p)- We do not, however, see strong 
evidence for the diffusive broadening of the wavefront predicted by ()14p . This could be because 
this equation applies only where the stationary density is small, a region that is a hard to access 
numerically. 

A key component of the scaling picture (section 4), seen explicitly within the mean- field 
dynamics (section 5), is the continuity of the intrinsic wave velocity across the transition. After 
unfolding the envelope (as in Fig. [7]) one can obtain this velocity over a range of activation 
rates r. In Fig. [8l this velocity is compared with the apparent velocity before unfolding the 
velocity (i.e., the data of Fig. [3|). The two velocities are clearly distinct in the subcritical regime 
[r < re ^ 3.3). Since the envelope decays exponentially in the subcritical phase, it is difficult 
to probe the late-time travelling-wave dynamics, and so the error bars on the modified wave 
velocity are necessarily large. However, the data suggest that the intrinsic velocity changes more 
slowly just below the critical point than above, in qualitative agreement with the mean-field 
predictions of equations p^ . lpOj) . 

Finally, we attempt to access the exponent X = ^11 ~ ^± by fitting the observed wave velocity 
in the supercritical regime to the form t;sup ~ Uc + A{r — Vc)^ suggested by the scaling picture of 
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Figure 8: The intrinsic velocity: Velocity vs. activation rate. We see a distinct difference 
between the apparent and intrinsic velocities. Note that the apparent velocities only 
seems continuous due to finite-time effects. 
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Figure 9: Velocity exponent x measured in the stochastic simulations. Shown is ln(w — Vc) 
against ln(r — re) with Vc — 1.795 to obtain a straight line in the regime r — >■ re- 
The solid line has the gradient x ~ 0.637 appropriate for Directed Percolation in one 
dimension. 



section 4 and confirmed within the mean- field regime by travelling wave analysis of section 5. It 
turns out that an estimate of x is rather sensitive to the values of Vc and Tc used in a straight-line 
fit to ln(z; — Vc) plotted as a function of ln(r — re). Taking re = 3.3055, we find the best fit (as 
quantified by the sum of square residuals) when Vc ~ 1.795. The corresponding plot is shown 
in Fig. El along with a line of gradient x = 0.637, which is the appropriate choice for DP in one 
dimension. We see that the simulation data have a gradient that is consistently higher than the 
predicted value. It may be that corrections to the leading-order behaviour remain significant in 
the region we have been able to access numerically. We remark that reasonable straight lines 
are obtained for values of Vc in the range 1.70 to 1.84, yielding estimates of x between 0.6 and 
0.8 suggesting that the numerical data are consistent with our scaling prediction for x- 

7 Conclusion 

In this work we have studied a variation of the contact process which includes an active boundary 
which drives the system and advection away from that boundary. As expected the system 
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exhibits a phase transition from a state with activity locahsed near the boundary to a state where 
a wave of activity emanates away from the boundary. We identify the DP critical exponents 
/3 and i^m by considering the behaviour of the density far from the boundary and the spatial 
decay length over which the density decays to that value. On the other hand, the velocity of 
the activity wave emanating from the boundary exhibits some perhaps unexpected behaviour: 
in the subcritical phase the apparent velocity is zero whereas in the supercritical phase the 
velocity jumps discontinuously to a non-zero value. We have explained this phenomenon by 
studying a mean field theory in which we show that an intrinsic velocity for a wave emanating 
from the boundary and described by ([7]), exists both below and above the transition. However 
in the subcritical regime the spatially decaying envelope for this wave means that the apparent 
velocity observed in the simulations is zero. This picture appears to hold well in simulations of 
the stochastic system 

The study raises several interesting questions. At the mean field level it would be interesting 
to put our analysis of the modified Fisher equation ([H]), in particular the subcritical case, on a 
more rigorous footing. A study of the possible crossover in the solution of Q from a nonlinear 
wave to a diffusive wave would be illuminating. 

The noisy version of the Fisher-KPP equation is known to describe contact processes p4j. 
It would of course be of great interest to further understand travelling wave solutions of the 
noisy version of the Fisher-KPP equation, in particular the velocity and width of the front. 

A Appendix: Solution of spatially discretized mean field equa- 
tion 



(21) 



The mean-field equation (jlSp governing the density in the DACP process is 

Pi = rpi^i{l- Pi)- Pi Vi > 1 

The steady-state solution {pi = 0) is 



(1 — r)r^ 
1 -r* 



P^ = V-IiTT (22) 



which yields the large i behaviour 



Pi ~ (1 - r)r^ for r < 1, (23) 

Pi = —— for r = 1, (24) 
I + 1 

Pi ^ ^"""^^ (l + r-^'+^A for r>l. (25) 



Thus this mean field theory is consistent with the profiles decaying a decay length which diverges 
as 1/1 lnr| ~ A~^ and the order parameter emerging as p ~ A where A = r — Vc with Vc = 1. 
Also at criticality we have a power law decay of the profile with exponent 5 = 1. All these 
exponents are consistent with the mean field DP exponents given in Table [TJ 
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